Lectures on Diffraction and Saturation of 



Nuclear Partons in DIS off Heavy Nuclei 1 

LP. Ivanov a , N.N. Nikolaev bc >, W. Schafer b \ B.G. Zakharov c ) 

& V.R. Zoller d ) 

a ) Institute of Mathematics, Novosibirsk, Russia 
b ) Institut f. Kernphysik, Forschungszentrum Jiilich, D-52425 Jiilich, Germany 
c > L.D. Landau Institute for Theoretical Physics, Chernogolovka, Russia 
d ) Institute for Theoretical and Experimental Physics, Moscow, Russia 



(N 
< 

o 



(N 
O 



E-mail: N.Nikolaev@fz-juelich. de 



Abstract 



The Lorentz contraction of ultrarelativistic nuclei entails a spatial overlap and fusion (recom- 
bination, saturation) of partons belonging to different nucleons at the same impact parameter. 
In these lectures we present a consistent description of the fusion of partons in terms of nuclear 
■ attenuation of color dipole states of the photon and collective Weizsacker- Williams (WW) gluon 



structure function of a nucleus. The point that all observables for DIS off nuclei are uniquely 
calculable in terms of the nuclear WW glue amounts to a new form of factorization in the sat- 
uration regime. 

We start with the theory of multichannel propagation of color dipoles in a nuclear medium 



including the color-singlet to color-octet to color-octet transitions. We show how the Glauber- 
Gribov formulas are recovered from the multichannel formalism. Then we derive the two-plateau 
momentum distribution of final state (FS) quarks produced in deep inelastic scattering (DIS) off 
nuclei in the saturation regime. The diffractive plateau which dominates for small p measures 
precisely the momentum distribution of quarks in the beam photon, the role of the nucleus is 
^ ■ simply to provide an opacity. The plateau for truly inelastic DIS exhibits a substantial nuclear 
broadening of the FS momentum distribution. The Weizsacker- Williams glue of a nucleus ex- 
hibits a substantial nuclear dilution, still soft initial state (IS) nuclear sea saturates because of 
the anti-collinear splitting of gluons into sea quarks. Then we comment on the signatures of 
saturation in exclusive diffractive DIS. 

A large body of these lectures is on the recent theory of jet-jet inclusive cross sections. We show 
that for hard dijets the decorrelation momentum is of the order of the nuclear saturation mo- 
mentum Qa- For minijets with the transverse momentum below the saturation scale we predict 
a complete disappearance of the azimuthal jet-jet correlation. We conclude with comment on 
a possible relevance of the decorrelation of jets to the experimental data from the STAR-RHIC 
Collaboration. 



1 The notes of lectures presented by N.N.N, at the XXXVI St. Petersburg Nuclear Physics Institute 
Winter School on Nuclear and Particle Physics & VIII St. Petersburg School on Theoretical Physics, 
St.Petersburg, Repino, February 25 - March 3, 2002 
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1 Introduction 



Within the standard QCD parton model the virtual photoabsorption cross section is 
proportional to the density of partons in the target. Inverting the argument, one would 
define the density of partons in a target A as 

F 2 (x,Q 2 ) = Y,e 2 f [q f (x,Q 2 ) + q f (x,Q 2 )} = j^—<r r A (1) 

f 47T Ct em 

In the so-called brick- wall or Breit frame, in which the photon has a zero energy and the 
(target) hadron or nucleus is ultrarelativistic, DIS amounts to backward reflection of the 
parton with the longitudinal momentum k z = xP z = \q z = |\/^F> see fig-1- 




Figure 1: DIS in the Breit frame 

When DIS is viewed in the laboratory frame, the hadronic properties of photons 
suggest P a nuclear shadowing by which the density of partons in a nucleus rises slower 
than oc A. Very heavy nuclei will be opaque to a part of the Fock states of the photon 
When viewed in the Breit frame, the Lorentz contraction of an ultrarelativistic nucleus 
entails a spatial overlap and fusion of partons at 

x< J x A = - (2) 

R A m N 

as indicated in fig. 2. 

This interpretation of nuclear opacity in terms of a fusion and saturation of nuclear 
partons has been introduced in 1975 pQ way before the QCD parton model. The review 
of early works on nuclear shadowing can be found in lectures by one of the authors at the 
1976 Winter School of PNPI see also The pQCD link between nuclear opacity 
and saturation has been considered in ref . [S] and by Mueller |B] , the pQCD discussion of 
fusion of nuclear gluons has been revived by McLerran et al. jjj- 

The common wisdom is that in DIS the final state (FS) interaction of the struck 
parton, which in the Breit frame moves in the direction opposite to the debris of the 
nucleon/nucleus, see fig. 2, can be neglected and the observed transverse and longitudinal 
momentum distribution of struck partons in the FS coincides with the initial state (IS) 
density of partons in the probed hadron. Based on the consistent treatment of intranuclear 



2 





-►x>x 




A No shadowing 



x « Shadowing 



Figure 2: The fusion of partons from overlapping parton clouds of two nucleons in a 
Lorentz- contracted nucleus 

distortions, we derive the two-plateau spectrum of FS quarks found in jH] . Following our 
early analysis of diffractive DIS off nucleons ED] and nuclei [TI], we introduce the 
collective Weizsacker- Williams gluon structure function of the nucleus, which uniquely 
defines all nuclear DIS observables. We find a substantial nuclear broadening of inclusive 
FS spectra and demonstrate that despite this broadening the FS sea parton density exactly 
equals the IS sea parton density calculated in terms of the WW glue of the nucleus as 
defined according to [TT]. We pay special attention to an important point that diffractive 
DIS in which the target nucleus does not break and is retained in the ground state, makes 
precisely 50 per cent of the total DIS events for heavy nuclei at small x j2] . We point out 
that the saturated diffractive plateau measures precisely the momentum distribution of 
(anti-)quarks in the qq Fock state of the photon. In contrast to DIS off nuclei, the fraction 
of DIS off free nucleons which is diffractive, and which also measures unitarity corrections 
to the two-gluon exchange approximation [T^IHH], is negligibly small [Hj, t}d ^ 6-10 %, 
and there is little room for genuine saturation effects even at HERA. We show how the 
anti-collinear splitting of WW gluons into sea quarks gives rise to nuclear saturation of 
the sea despite a substantial nuclear dilution of the WW glue. 

After establishing the properties of the single-jet inclusive cross section, we turn our 
attention to the recent theory of jet-jet (de) correlations [T^]. After a hard scattering par- 
tons with large transverse momentum fragment into a high-energy jet of particles. In deep 
inelastic scattering at small x, to the collinear approximation, the photon-gluon fusion 
1*9 ~ > QQi often referred to as an interaction of the unresolved photon, gives rise to the 
back-to-back jets in the current or photon fragmentation region. The experimental signa- 
ture of the unresolved photon interaction is that the so-called lightcone plus- components 
of the jet momenta sum up to the lightcone plus-component of the photon's momentum. 
The allowance for the transverse momentum of gluons leads to the disparity of the mo- 
menta and azimuthal decorrelation of the quark and antiquark jets, which within the 
/c^-factorization can be quantified in terms of the so-called unintegrated gluon structure 
function of the target (see ^H] and references therein). 
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In view of the substantial nuclear broadening of the unintegrated gluon SF of nuclei it 
is natural to expect stronger azimuthal decorrelation of jets produced in DIS off nuclei. In 
these lectures we report a consistent multiple-scattering theory of jet-jet decorrelation in 
the genuinely inelastic DIS followed by color excitation of the target nucleus. For heavy 
nuclei of equal importance is coherent diffractive DIS in which the target nucleus does 
not break and is retained in the ground state. In these coherent diffractive events quark 
and antiquark jets are produced exactly back-to-back with negligibly small transverse 
decorrelation momentum 

|A| = |p+ + p_| < 1/R A ~ mJA 1 ^ , (3) 

where Ra is the radius of the target nucleus. For hard jets diffractive attenuation effects 
are weak and we report a compact formula for the broadening of azimuthal correlations 
between the quark and antiquark jets. In this case the decorrelation (acoplanarity, out- 
of-plane) momentum is given by the nuclear saturation momentum Qa- Then we prove a 
disappearance of the jet-jet correlation for minijets with the transverse momentum below 
the saturation scale Qa- In the Conclusions we summarize our principal findings and 
speculate on their relevance to the recent observation of the dissolution of the away jets 
in central nuclear collisions at RHIC |17j . 



2 Quark and antiquark jets in DIS off free nucleons: 
single particle spectrum and jet-jet decorrelation 

We recall briefly the color dipole formulation of DIS [21 El 13 UH U21 UH| and illustrate our 
ideas on an example of jet-jet decorrelation in DIS off free nucleons which at moderately 
small x is dominated by interactions of qq states of the photon. We evaluate the jet 
production cross section at the parton level. The total cross section for interaction of the 
color dipole r with the target nucleon equals 

cr(r) = a s (r)a J d 2 nf(hi) [1 - exp(«Kr)] 

= 7}®s(r)a J d 2 Kf(K) [1 - exp(iKr)] • [1 - exp(-iicr)] , (4) 

where /(k) is related to the unintegrated glue of the target nucleon by 

/( } ~N c a ^ dlog* 2 [b) 

and is normalized as 

>*/(«) = !. (6) 



For DIS off a free nucleon target, see figs. 3a-3d, the total photoabsorption cross 
section equals 

a N = [ d 2 rdz\ty(z,r)\ 2 a(r) . (7) 
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Figure 3: The pQCD diagrams for inclusive (a-d) and diffractive (e,f) DIS off protons 
and nuclei (g-k). Diagrams (a-d) show the unitarity cuts with color excitation of the 
target nucleon, (g) - a generic multiple scattering diagram for Compton scattering off 
nucleus, (h) - the unitarity cut for a coherent diffractive DIS, (i) - the unitarity cut for 
quasielastic diffractive DIS with excitation of the nucleus A* , (j,k) - the unitarity cuts for 
truly inelastic DIS with single and multiple color excitation of nucleons of the nucleus. 
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Upon the relevant Fourier transformations one finds the momentum spectrum of the final 
state (FS) quark prior the hadronization, 

da N a a s (p 2 ) f 2 f r \ \/ *\ \ i *i \|2 / c \ 

d^d~ z = y --fry 1 s{) 1(7 |p) " (7 |p " K)l {) 

where p is the transverse momentum of the quark and z + = z and z_ = 1 — z are the 
fractions of photon's lightcone momentum carried by the quark and antiquark, respec- 
tively. The variables z± for the observed jets add up to unity which is the signature of 
the unresolved photon interaction. 

For the reference, we cite here the relevant overlap of the momentum space wave 
functions of the photon, summed over the helicities of the FS quark and antiquark [Sj. 
For transverse photons 
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and for the longitudinal photons 
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where 
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e 2 = z(l - z)Q 2 + mj 

and my is the mass of the quark qf with the charge e/. Now, notice that the transverse 
momentum of the gluon is precisely the decorrelation momentum A, so that in the still 
further differential from 

>2\ 



da 



N 



a a s {p 



dzd 2 p + d 2 A 2 (2tt) 2 
A useful small-A expansion for light flavors is 

da N N c e 2 f a em a a s (p 2 ) 
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(12) 



dzd 2 p + d 2 A (2tt) 2 
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(13) 



where and Aj_ are the components of the decorrelation parallel, and transverse, to p + . 
Only the latter, acoplanarity component, contributes to the azimuthal jet-jet decorrela- 
tion, the former is only measurable if full experimental reconstruction of jets is possible. 
This is the leading order result from the k± factorization, for the applications to DIS off 
free nucleons see ^B] and references therein. Below we shall extend (jl2j) to jet-jet inclusive 
cross section for interactions of unresolved photons with nuclear targets. 

3 Propagation of color dipoles in nuclear medium 

We focus on DIS at x < xa = ^/Ra^n "C 1 which is dominated by interactions of 
qq states of the photon. Extension to interactions of higher Fock states of the photon 
is straightforward and will give rise to the evolution effects, which we do not consider 
here. For x < xa the propagation of the qq pair inside nucleus can be treated in the 
straight-path approximation. 

We work in the conventional approximation of two t-channel gluons in DIS off free 
nucleons. The relevant unitarity cuts of the forward Compton scattering amplitude are 
shown in fig. 3a-3d and describe the transition from the color-neutral qq dipole to the 
color-octet qq pair. (Really, for SU(N C ) we have fundamental & adjoint multiplets instead 
of the triplet & octet familiar for N c = 3, our continuous reference to the latter should not 
cause any confusion.) The two-gluon exchange approximation amounts to neglecting small 
unitarity corrections to DIS off free nucleons [T2*l I13j. which are quantified by diffractive 
DIS described by higher order diagrams of fig. 3e,3f. This approximation is justified by a 
small fraction of diffractive DIS, f|£) <C 1 [Hj. The unitarity cuts of the nuclear Compton 
scattering amplitude which correspond to the genuine inelastic DIS with color excitation 
of the nucleus are shown in figs. 3j,3k. Here we notice that in the diagram 3k the qq pair 
propagates in the both color-octet and color-singlet states. 

Let b + and b_ be the impact parameters of the quark and antiquark, respectively, 
and 5'yi(b + ,b_) be the S-matrix for interaction of the qq pair with the nucleus. We are 
interested in the inclusive cross section when we sum over all excitations of the target 
nucleus when one or several nucleons have been color excited. A convenient way to sum 
such cross sections is offered by the closure relation. Regarding the color states c^ m of 
the qkq m pair, we sum over all octet and singlet states. Then, the the 2-body inclusive 
spectrum is calculated in terms of the 2-body density matrix as 



Here \1/ is the gg-Fock state wave function of the virtual photon, and we suppressed its 
dependence on z±. In the integrand of (|14j) we subtracted the diffractive component of the 





x exp[— zp + (b + — b' + ) — ip_(b_ — b'_)] 
x^*(b' + -b'_)^(b + -b_) 

x {e EU; A\s* A Q>+> b '-)l^*; c *-> ( c *™; a*\SaQ>+, b 



(14) 
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final state. Notice, that the four straight-path trajectories b±, b' ± enter the calculation of 
the full fledged 2-body density matrix and Sa and S A describe the propagation of two qq 
pairs inside a nucleus. 

Upon the application of closure to sum over nuclear final states A* the integrand of 
()14j) can be considered as an intranuclear evolution operator for the 2-body density matrix 
(for the related discussion see ref. |20j) 

EE(^l{(l|^(bV,b'_)|c fcm )}|A*)(A*|{( Cfem |^(b + ,b_)|l)}|A) = 

A* km 

= (A\ [x)<l|5 , I(b;,b , _)|c fcro )(c Jfcm |5' A (b + ,b_)|l)] \A) (15) 

k km ) 

Let the QCD eikonal for the quark-nucleon and antiquark-nucleon one-gluon exchange 
interaction be T?A(b) and T°A(b), where Ti and T" are the SU(N C ) generators for 
the quark and antiquark states, respectively. The vertex V a for excitation of the nucleon 
g a N — > N* into color octet state is so normalized that after application of closure the 
vertex g a g b NN in the diagrams of fig. 3a-3d is 5 a b. Then, to the two-gluon exchange 
approximation, the S-matrix of the (gg)-nucleon interaction equals 

S N (b+, b_) = 1 + z[T»A(b+) + r«A(b_)]F - i[T^A(b + ) + T*A(b_)] 2 (16) 

and the color-dipole cross section for the qq dipole is 

N 2 — 1 r 

a(b + - b_) = H^- J d 2 b + [A(b + ) - A(b_)] 2 . (17) 
The nuclear S'-matrix of the straight-path approximation is given by 

A 

S A {b+, b_) = J] S N (h+ - b jt b - bj) (18) 

where bj are the transverse coordinates of nucleons in a nucleus and the ordering along 
the longitudinal path is understood. We evaluate the nuclear expectation value in (|15|) in 
the standard dilute gas approximation. To the two-gluon exchange approximation, only 
the terms quadratic in A(b,) must be kept in the evaluation of the single- nucleon matrix 
element 

(A^(b; - b if b'_ - b 3 )S N (b + - bj, b - b^Nj) 

per each and every nucleon which enters the calculation of S\Sa- Following the technique 
developed in [2UI22] we can reduce the calculation of the evolution operator for the 2-body 
density matrix (|T5)l to the evaluation of the S'-matrix S 4 ^(b + , b_, b' + , b'_) for a scattering 
of a fictitious 4-parton state composed of the two quark-antiquark pairs in the overall 
color-singlet state. Because (T")* = — Tf, within the two-gluon exchange approximation 
the quarks entering the complex-conjugate S A in (|T5)l can be viewed as antiquarks, so 
that 

({ y Z(l\S* A (b' + ,b'_ )|c fcm ><c fcm |S A (b +) b_)|l> 

I km 

= ^2 S k i8 mj (c km Cji\S4A(b' + , b'_,b + ,b-)\ll) (19) 

kmjl 




s 



where S^ib^, b'_, b + , b ) is an S-matrix for the propagation of the two quark-antiquark 
pairs in the overall singlet state. While the first qq pair is formed by the initial quark 
q and antiquark q at impact parameters b + and b_, respectively, in the second pair q'q' 
the quark q' propagates at an impact parameter b'_ and the antiquark q' at an impact 
parameter b^_, as indicated in fig. 4. In the initial state the both quark-antiquark pairs 
are in color-singlet states: \in) = |11). 
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Figure 4: The pQCD diagrams for the matrix of color dipole cross section for the J^-body 
{qq)(q'q') state. The set J^i-J^l shows only half of the diagrams for scattering with rotation 
of the color state of dipoles. 



Let us introduce the normalized singlet-singlet and octet-octet states 

111) = ^(W)(§V) 



{qT a q){q'T a q" 



N*-l 



(20) 
(21) 



9 



where N c is the number of colors and T a are the generators of SU(N C ). Making use of 
the color Fiertz identity, 

dkmdji = ±6fS? + 2 J2(T a )t(T a )T ( 22 ) 

iV c a 

the sum over color states of the produced quark-antiquark pair can be represented as 
|S 4A (b' + ,b'_,b + ,b_)|ll> = (ll|^(b' + ,b / _,b + ,b_)|ll) 

km 

+ v /iV c 2 -l(88|5 4A (b / +! b'_ ! b + ,b_)|ll) (23) 

If cr 4 (b' + , b'_, b + , b_) is the color-dipole cross section operator for the 4-body state, 
then the standard evaluation of the nuclear expectation value for a dilute gas nucleus 
gives 

S 4A {b' + , b'_, b +) b_) = exp{-ia 4 (b' + , b'_, b + , b_)T(b)} (24) 

where 

T(b) = J dzn A (z,b) (25) 

is the optical thickness of a nucleus at an impact parameter b rs h±,b' ± . In writing down 
()19|) we employed the standard approximation of neglecting the size of color dipoles com- 
pared to a radius of heavy nucleus. The single-nucleon S'-matrix contains transitions 
from the color-singlet to the both color-singlet and color-octet qq pairs. However, only the 
color-singlet operators contribute to (Nj\S^(h' + — bj, b'_ — bj),Sjv(b+ — bj, b_ — bj)\Nj), 
and the matrix <r 4 (b+, b'_, b + , b_) only includes transitions between the |11) and 1 88) 



(26) 

-K) 

- K)} (27) 

b'J] 
b' + )] 

b'_)] (28) 



color-singlet 4-parton states, the 1 18) states are not allowed. 
Performing the relevant color algebra, we find 

(11|ct 4 |11) = a(b + - b_) + a(b' - b' 



(11|<7 4 I88)= . [ <r(b_ - b'J + a(b^ 

— cr(b + — b'__) — cr(b_ 



N 2 — 2 

(881^88) = -f—i[a(b + - b' + ) + a(b_ 

c 

' >(b+ - b_) + a(b' 



N? - 1 
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The term in (|14j) . which subtracts the contribution from processes without color excitation 
of the target nucleus, equals 



(A\ {<l|Sl(b' +) b'_)|l>} \A)(A\ {(l|fiU(b +I b_)|l» \A) 



exp{-- [a(b+ 



a(b' + -W_)] T(b)} 



(29) 



Let Ei j2 be the two eigenvalues of the operator 04. Then it is convenient to use the 
Sylvester expansion 



exp{-ia 4 T(b)} = 
eXp{ -^ SlT(b)} |r^ 



^ + exp{-^ 2 T(b)}|^| 



(30) 



Evidently, in the general case the eigen-cross sections are fairly complex functions of the 
six possible dipoles formed by four partons. One may wonder how do the standard color- 
dipole Glauber-Gribov formulas |2| for nuclear DIS cross sections emerge from ([30)1? 
Hence we turn first to the simpler cases of the total inelastic and single-jet inclusive cross 
sections. An application to of the Sylvester expansion gives 



((11| + v/iV c 2 -l(88|) exp {-^(b)} |11) 

-exp{-i[a(r) + ( x(r')]T(b)} 

= exp {-^S 2 T(b)} - exp |-± |a(r) + a(r')} T(b) 



, (ll|a 4 |ll>-E 2 
+- — — ^ exp 



Ei — E2 



exp 



--E 2 T(b) 



- l(ll|o- 4 |88) 
Ei — E 2 



exp 



;SiT(b) 



exp 



-E 2 T(b) 



(31) 



4 Non- Abelian aspects of propagation of color dipoles 
in nuclear medium and Glauber-Gribov formalism 

In order to get an insight into the impact of propagation of the color-octet qq pairs inside 
a nucleus let us consider first the total inelastic cross section obtained from ()14|) upon 
the integration over the transverse momenta p± of the quark and antiquark. Such an 
integration amounts to putting b + = b' + and b_ = b'_. Then we are left with the system 
of two color dipoles of the same size r = b + — b = r' = b' + — b'_, and the matrix of the 
4-body cross sections has the eigenvalues 



Ei = 0. 



(32) 



N? - 1 



a(r) 



(33) 
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with the eigenstates 



IA) = ^(|H) + V /iV c-l|88)), (34) 



|/ 2 > = — (\/iV c 2 -l|ll>-|88>). (35) 

The eigen-cross section (}3*3*j) differs from cr(r) for the color-singlet qq pair by the nontrivial 
color factor N 2 /(N 2 — 1). The existence of the non-attenuating 4-quark state with Si = 
is quite obvious and corresponds to an overlap of two dipoles of the same size with 
neutralization of color charges. 

Now notice, that the final state which enters the calculation of the genuine inelastic 
DIS off a nucleus, see eq. (J23|) . is precisely the eigenstate Then, even without 

invoking the Sylvester expansion ()3(Jj). ()31j). the straightforward result for the inelastic 
cross section is 



x 



J d 2 rdz\V(z,r)\ 2 



-ia 4 T(b) 



|11) -exp [-a(r)T(b)] 



= J d 2 b(l*\ {exp [™£iT(b)] - exp (-a(r)T(b)]} | 7 *) 

= J d 2 b( 7 *| {1 - exp [-a(r)T(b)]} | 7 *> (36) 

what is precisely the Glauber-Gribov formula [2] in which no trace of a propagation inside 
a nucleus of the eigenstate (fSHj) with the eigen-cross section is left. 

When the final state qq is produced in the color-singlet state, the net flow of color 
between the qq pair and color-excited debris of the target nucleus is zero. Which suggests 
that a rapidity gap may survive upon the hadronization, although whether a rapidity 
gap in genuine inelastic events with color-singlet qq production is stable against higher 
order correction or not remains an interesting open issue. Here we proceed under the 
assumption that the rapidity gap survives and the production of the color-singlet and 
color-octet qq pairs can be separated. Although the debris of the target nucleus have 
a zero net color charge, the color-excited debris of nucleons are spatially separated by 
a distance of the order of the nuclear radius. The formation of color strings between 
the color centers would lead to a total excitation energy of the order of 1 GeV times 
A 1 / 3 , so that such a rapidity-gap events would look like a double diffraction with multiple 
production of mesons in the nucleus fragmentation region (for the theoretical discussion of 
conventional mechanisms of diffraction excitation of nuclei in proton-nucleus collisions see 
|23| . the experimental observation has been reported in [21] )■ As such it is distinguishable 
from quasielastic diffractive DIS followed by excitation and breakup of the target nucleus 
without production of secondary particles. 

Making use of the Sylvester expansion (|3*U J) -([3~T |) . and the eigenstates (|53 jl .(|35 |) . one 
readily obtains 

a in (A*(qq)i) = 



12 



exp [-<7(r)T(b)]) 



iV 2 



N* 



[1 — exp 



1 



£ 2 T(b) 



It*) 



and 



a in (A*(gg) 8 ) 



iV 2 



d 2 b(Y 



exp 



--E 2 T(b) 



17 ) • 



(37) 



(3? 



In the latter case the nontrivial eigen-cross section X 2 enters explicitly. 

Several features of the result (J37|) are noteworthy. First, the color neutralization of 
the qq pair after the first inelastic interaction requires at least one secondary inelastic 
interaction. Indeed, an expansion of the integrand starts with the term quadratic in the 
optical thickness: 



(1-exp [— <x(r)T(b)]) 

q 2 (r)T 2 (b ) 

= 2(iV 6 2 - 1) 



1 



exp 



-\^T{h) 



+ ... 



(39) 



Second, in the large- iV c limit the color octet state tends to oscillate in color remaining in 
the octet state. This is clearly seen from (J39|) . Third, in the limit of an opaque nucleus 



a in (A*(qq)i) = ^ J d 2 h(l*\ {1 " exp [-<r(r)T(b)]} | 7 *> 



(J; 



(40) 



and remains a constant fraction of DIS in contrast to the quasielastic diffractive DIS 
or inelastic diffractive excitation of a nucleus, the cross sections of which vanish for an 
opaque nucleus j2 • 

The result ()38|) shows how the non-Abelian nature of the propagation of color dipoles 
inside a nucleus which manifests itself via the eigen-cross section (}3*3*j) distinct from cr(r) 
reveals itself in the total cross section of inelastic DIS with excitation of dijets in the 
color-octet state. 

An analysis of the single-parton, alias the single-jet, inclusive cross section is quite 
similar. In this case we integrate over the momentum p_ of the antiquark jet, so that 
b'_ = b_. The matrix o" 4 has the eigenvalues 



Si 



a r 



(41) 



N tj Hr) + <x(r')] - ^rr^r - (42) 



AT2 

with exactly the same eigenstates and |/ 2 ) as given by eqs. (|33 jl .(J35 |) . Again, the 
cross section of genuine inelastic DIS corresponds to projection onto the eigenstate 
so that EI 



da i n 



d 2 hd 2 pdz (2tt 



^— / tZ 2 r'd 2 rexp[ip(r'-r)]**(r')*(r) 

Tl) 2 J 
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x {expf-i^TCb)] - exp[-~[a(r) + <x(r')]T(b)]} 
= — -- / d 2 r'd 2 rexp[ip(r' - r)]**(r')#(r) 

(27TJ 2 J 

x {exp[-^(r - r')T(b)] - exp[-±[<r(r) + a(r')]T(b)]} . (43) 

5 The Pomeron-Splitting Mechanism for Diffractive 
Hard Dijets and Weizsacker- Williams glue of nuclei 

At QCD parton level diffraction dissociation of photons and hadrons is modeled by exci- 
tation of qq, qqg, ... Fock states which are lifted on their mass shell through the t-channel 
exchange of a QCD Pomeron with the target hadron. The color singlet two-gluon struc- 
ture of the Pomeron gives rise to the two distinct forward qq dijet production subprocesses: 
The first one, of fig. 3e (also fig. 5a), is a counterpart of the classic Landau-Pomeranchuk 
|2H1 EEl EZ| mechanism of diffraction dissociation of deuterons into the proton-neutron 
continuum and can be dubbed the splitting of the beam particle into the dijet, because 
the transverse momentum k of jets comes from the intrinsic transverse momentum of 
quarks and antiquarks in the beam particle. Specific of QCD is the mechanism of fig. 3f 
(also fig. 5b) where jets receive a transverse momentum from gluons in the Pomeron. In 
an extension of the early work j^j, in [lOj it was shown that the second mechanism dom- 
inates at a sufficiently large k and in this regime diffractive amplitudes are proportional 
to the differential (unintegrated) gluon structure function jF(x,k 2 ) = dG(x, k 2 )/<91ogk 2 . 
Correspondingly, this mechanism has been dubbed splitting of the Pomeron into dijets. 
In diffractive DIS the Pomeron splitting dominates at k ^> Q, whereas the somewhat 
modified Landau et al. mechanism dominates at k < Q. 

Motivated by the recent data from the E791 Fermilab experiment [2E]> nere we discuss 
peculiarities of the Pomeron splitting mechanism on an example of diffractive excitation 
of pions into dijets on free nucleon and nuclear targets. In a slight adaption of the results 
of ED] one should replace the pointlike j*qq vertex eA^^^ by the non-pointlike nqq 
vertex iT{M 2 )^!'^^> . Here M 2 = (k 2 + m 2 )/z(l — z) is the invariant mass squared of the 
dijet system, the pion momentum is shared by the jets in the partitioning z, 1 — z. The 
explicit form of the spinor vertex reads 

^ A (k) 75 ^A(-k) = A [m f 5 x _- x - V2k ■ e_Ax] , (44) 
' z(l — z) 



where is the quark mass, A and A are the quark and antiquark helicities and e x = 
— (\e x + ie y )/y/2. The two helicity amplitudes &o(z, k, A) for A + A = and &i(z, k, A) 
for A + A = ±1, can be cast in the form (we concentrate on the forward limit, A = 0) 



$o(z, k) = a s (k 2 )a 



J d 2 nmfijj 7T (z,'k)f(K) — J d 2 K,mf^jj 7T (z, «)/(k 



(45) 
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and $1 is obtained from the substitution rrifip(z, k) — > \£ip(z, k). The radial wave function 
of the qq Fock state of the pion is related to the nqq vertex function as 

^ k) = ^ Z (l-z) (M>- ml) (46) 
and is normalized to the ir —>■ fiu decay constant F n = 131 MeV through 

F w = J d 2 'kdzmf^p 7V (z,'k) = J dz^(z) , (47) 

where ifi^z) is the pion distribution amplitude Finally, our normalization of helicity 
amplitudes is such, that the differential cross section of forward dijet production equals 



da 



D 



dzdk 2 dA 2 



A=o 24 



^{l$o| 2 +|*i| 2 }. (48) 



We turn to the discussion of the asymptotics for large jet momenta k. The first term 
in eq. (|45j) comes from the Landau et al. pion splitting mechanism of fig. 5a, whereas the 
second one is the contribution from the Pomeron splitting of fig. 5b. Because ip^ <g) / is a 
convolution of nonoscillating functions, it is necessarily a broader function than ip w (z, k) 
and thus will take over if only k is large enough. For the quantitative estimate, we remind 
the reader, that for x ~ 10~ 2 relevant to the kinematics of E791, the large-k behavior of 
/(k) is well described by the inverse power law /(k) oc k~ 2S with an exponent 5 ~ 2.15 
|3Uj . Clearly, /(k) decreases much slower than ifj(z,k), and hence the asymptotics of the 
convolution integral is controlled by the asymptotics of /(k): 

J d 2 Km f ip 7T (z, K)f(k - k) w f(k) J d 2 Km f 'ilj 7r (z 1 k) = f(k) ^(z) F n , (49) 

which shows that in this regime the dijet momentum comes entirely from the momentum 
of gluons in the Pomeron. Furthermore, in the same regime diffraction into dijets probes 
the pion distribution amplitude ^(z). 

Now we notice that in the color dipole representation 

$oO, k) = J d 2 r e~ lkr cr(r) m f ^ T (z, r) , (50) 

and the nuclear amplitude is readily obtained [S] by substituting in eq.(f5TH) 

<7(r) ^ a A {v) = 2 J d 2 b |l - exp[-^(r)T A (b)]} (51) 

Typical nuclear double scattering diagrams of figs. 5c-5e can conveniently be classified as 
nuclear shadowing of the pion splitting (fig. 5c), nuclear shadowing of single Pomeron 
splitting (fig. 5d) and double Pomeron splitting (fig. 5e) contributions. The j Pomeron 
splitting is due to diagrams, in which j of the Pomerons exchanged between the color 
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p - (i-z) A 

x i\ C x 2 P* +Z A 
k - A/2< < k + A/2 





c) 



d) 



e) 



Figure 5: Sample Feynman diagrams for diffractive dijet excitation in tcN collisions [di- 
agrams 5a), 5b)] and typical rescattering corrections to the nuclear coherent amplitude 
[diagrams 5c),5d),5e)}. 

dipole and the nucleus couple with one gluon to the quark and with one gluon to the 
antiquark and involve the j-fold convolutions of the unintegrated gluon distribution 



3 

E 

i=l 1=1 

Now we can invoke the NSS representation for the nuclear attenuation factor flj |S] 
1 



(52) 



exp 



-a(s)T(b) 



exp \-v A {h)} exp 



vaQd) J d 2 nf{n) exp(zAts) 
J d 2 K$(v A (b), k) exp(ms) (53) 
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where 



$(u A (b), k) = exp(-i/ A (b))/(°)(K) + (farwQ>, k) 



(54) 



and it is understood that 



vaO>) = -a s (r)a T(b) 



/<°>(k) = <$(«)• 



(55) 



(56) 



Then the multiple-Pomeron splitting expansion for nuclear diffractive amplitude takes the 
form 



$^ ) (^,k,A) = 2m / [d 2 be- thA 



i>i 



>ww 



(k-K) 



(57) 



A comparison with (J4~5j) suggests that l yw( K ) can be identified with the unintegrated 
nuclear Weizsacker- Williams glue per unit area in the impact parameter plane [8 J . In (|54|) 
we kept an explicit dependence on the optical thickness of the nucleus. 



6 Nuclear dilution and broadening of the unintegrated 
Weizsacker- Williams glue of nuclei 

On the one hand, according to ^IJ |H] the multiple convolutions f^(n) have a meaning 
of the collective unintegrated gluon structure function of j nucleons at the same impact 
parameter the Weizsacker- Williams gluon fields of which overlap spatially in the Lorentz 
contracted nucleus. On the other hand, these convolutions of can also be viewed as a 
random walk in which /(/«) describes the single walk distribution. To the lowest order in 
pQCD, the large k 2 behavior is 



/(«) oc 



2\2 



(58) 



The QCD evolution effects enhance /(/«) at large k 2 , the smaller is x the stronger is 
an enhancement. Because /(k) decreases very slowly, we encounter a manifestly non- 
Gaussian random walk. For instance, as was argued in JT], a j-fold walk to large k 2 is 
realized by one large walk, n\ ~ k 2 , accompanied by (j — 1) small walks. We simply cite 
here the main result ^T] 



/W(#c) =j-f(n) 



4tt 2 (j - 1)Y 
N c a K 2 



■ G(n 2 ) 



(59) 
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where G(k 2 ) is the conventional integrated gluon structure function and 7 is an exponent 
of the large-K 2 asymptotics of the differential glue in the proton, 



/(«) °c 



2Yy 



(60) 



Then the hard tail of unintegrated nuclear glue per bound nucleon can be calculated 
parameter free: 



fww(*>) = f( K ) 



1 + n at 2 G ( K ' 



(61) 



In the hard regime the differential nuclear glue is not shadowed, furthermore, because of 
the manifestly positive-valued and model-independent nuclear higher twist correction it 
exhibits nuclear antishadowing property [TT] . 

Now we present the arguments in favor of the soft-K behavior 



Q 2 i 



1 n 2 



with the width 



7t{k 2 + Q 2 ) 2 QjQ 



Q) ~ jQl 



In an evolution of f^\n) with j at moderate rc 2 , 

/°' +1 V)= /rf 2 k/(k 2 )/^((^-k) 2 ) 



(62) 



(63) 



(64) 



the function /(k 2 ) is a steep one compared to a smooth and broad /^((k — k) 2 ), so that 
we can expand 



K 



d* 2 



■[k 2 -2«k] + - 



lc^'V; 



/ {J V) + 
/ (i V) + 



d* 2 



dn 2 



dn 2 



(dK 2 ) 2 

k 2 



2 (dK 2 ) 2 

k 2 



4(/dc) 2 



(65) 



Here ==>■ indicates an azimuthal averaging. Our expansion was valid for k 2 ^Q], hence 



/(, +1 V) = ^V) + 



d 

di^ 2 



Kj~ 



dn 2 



d 2 k/(k 2 )k 2 



/^V) + <?0') 



rf 
d^ 2 



d^ 2 



(66) 
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where 

Arc 



/V ■ 4-7T 
3 d 2 kk 2 /(k 2 ) = j^G^Q)) (67) 



N c a 

is a smooth function of j. 

For small k 2 the recurrence relation amounts to 



which, for large j, reduces to a differential equation 



$ x(o) 

with the solution 



g(j) (69) 



,2 . X'(0)) fi , „• , -nX'(O) 



« = -^f y W) - -^o-)^ - j<& (70) 

The expansion (J55j) holds up to the terms oc k 2 and its differentiation at n 2 = gives a 
similar constraint on the j-dependence of Q 2 . The soft parameters Ql and o"o are related 
to the integrated glue of the proton in the soft region, 

2 271 " 2 

Qo a o ~ -j-j-G so ft , G so /t ~ 1 • (71) 

We conclude by the observation that when extended to n 2 > Q 2 , the parameterization 
(JE2J), (jHHl) has the behavior jQl/(K 2 ) 2 which nicely matches the j-dependence of the 
leading twist term in the hard asymptotics (jfiljl . 

An approximation ()62|). (|63j) is corroborated by the numerical results for the j-dependence 
of convolutions calculated for the unintegrated gluon structure function of the proton from 
ref. jSD], see fig. 6. 

For heavy nuclei the dominant contribution to 4>ww{k) comes from j « u A (b), and 
the gross features of the WW nuclear glue in the soft region are well reproduced by 

^ K) ™\ {K*%tf > (72) 

where the saturation scale Q\ = va(^)Qq oc A 1 / 3 . 
Notice a strong nuclear dilution of soft WW glue, 

4)ww{k) oc l/Q\ oc A' 1 ' 3 (73) 

which must be contrasted to the A-dependence of hard glue, cf. ()59|). (JfiT]) . 

<P ww (k) = u A (b)fww(K) oc A 1/3 x (1 + A 1/3 x (#T) + ...) . (74) 

Take the platinum target, A = 192. The numerical estimates show that for (qq) color 
dipoles in the average DIS on this target v& ~ 4 and Q 2 A ~ 0.8GeV 2 . For the qqg Fock 
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Figure 6: The dilution for soft momenta and broadening for hard momenta of the multiple 
convolutions /^(k) • 



states of the photon, which behave predominantly like the dipole made of the two octet 
color charges, the saturation scale is larger by the factor Ca/Cf = 9/4, 

QIa = ~ptQ\a ■ (75) 
The numerical estimates for the the average DIS on the Pt target give Ql A ~ 2.2 GeV 2 . 

7 Back to the hard diffraction of pions: the interpre- 
tation of the E791 data 

The detailed numerical analysis for the kinematics of the E791 experiment is found in 
In fig. 7 we only show the numerical result of the k dependence of the dijet cross 
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section with the (unnormalized) data from E791. We note in passing that the region of jet 
momenta k < 1.5 GeV is contaminated by diffractive excitation of heavy mesons ai,7r', 
etc, and in this region the use of plane wave parton model formulas is not warranted. Our 
calculations for smaller k only serve to give an idea on how the pion splitting dominates 
at small k and how the pomeron splitting mechanism takes over beyond the dip in $q 
and The latter is a higher twist term. In the pomeron splitting dominance region of 
k > 1.5 GeV we find good agreement with the experimentally observed k-dependence. We 
wish to warn the reader however, that for the kinematics of E791 the dijet cross section 
receives a huge contribution from the antishadowing nuclear higher twist correction [TT] . 




k [GeV] 



Figure 7: The E791 data I28f for the differential diffractive dijet cross section da/dk for 
the 196 Pt target with the theoretical calculations. The data are not normalized. The dash- 
dotted line shows the contribution of the helicity amplitude & A \ the dashed line is the 
contribution from <&[ A \ The solid line is the total result. 

The effect of antishadowing nuclear higher twist correction is especially obvious in 
fig. 8, where we show the exponent a of the A-dependence o ~ A a . In the impulse 
approximation 

° « 4- ~ a4/3 ( 7e ) 

K A 

If the exponent a is defined from a comparison of the cross sections for the carbon and 
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platinum targets, then the impulse approximation gives a ~ 1.44. For hard dijets our 
antishadowing effect predicts a substantial excess of a over the impulse approximation 
value. 
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Figure 8: The E791 data [28] for the exponent of the A dependence of the jet cross section 
vs. the theoretical calculations shown by thge solid curve. 



The NSS dominance of the pomeron-splitting contribution for hard dijets has fully been 
confirmed by the collinear factorization NLO order analysis of Chernyak et al. [HI] and 
Braun et al. [12] • The NLO correction to the NSS amplitude is found to be proportional 
to the asymptotic distribution amplitude and numerically quite substantial, so that the 
experimental data by E791 [2H] can not distinguish between the asymptotic and double- 
humped distribution amplitudes. According to NSS [TT] realistic model distributions 
do not differ much from the asymptotic one, though. To our view, the only caveat in 
the interpretation of the NLO results is that the issue of partial reabsorption of these 
corrections into the evolution/renormalization of the pion distribution amplitude still 
needs further scrutiny. Anyway, there emerges a consistent pattern of diffraction of pions 
into hard dijets and in view of these findings the claims by Frankfurt et al. [HH] that 
the diffractive amplitude is proportional to the integrated gluon structure function of the 
target must be regarded null and void. (Hopefully, some day the E791 collaboration shall 
report on the interpretation of their results within the correct formalism.) 
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8 Nuclear partons in the saturation regime 



Here we apply the results of sections 4 and 5 to the calculation of FS single-parton (jet) 
spectrum for nuclear targets. On the one hand, making use of the NSS representation, 
the total nuclear photoabsorption cross section (|51jl can be cast in the form 

a A = Jd 2 bJdzJ-0^J d 2 K<p ww (K) |< 7 *|p> - ( 7 *|p - k)\ 2 (77) 

which has a profound semblance to (JHJ) and one is tempted to take the differential form 
of (|77jl as a definition of the IS sea quark density in a nucleus (hereafter we consider one 
flavour and take e 2 the factor out of the photon's wave function and cross sections): 

dqis 1 Q 2 da A ^ 



d 2 hd 2 p 2 47r 2 a em d 2 hd 2 p 

In terms of the WW nuclear glue, all intranuclear multiple-scattering diagrams of fig. 3g 
sum up to precisely the same four diagrams fig. 3a-3d as in DIS off free nucleons. Further- 
more, one can argue that the small-x evolution of the so-defined IS nuclear sea is similar to 
that for a free nucleon sea. Although p emerges here just as a formal Fourier parameter, 
we shall demonstrate that it can be identified with the momentum of the observed final 
state antiquark. 

On the other hand, making use of the NSS representation, after some algebra one finds 
exp[ _I a(r _ r ') T(b )] - exp[-±[a(r) + a(v')}T(b)} = 
J d 2 K(p ww (K,){(exp[iK,(r — r')] — 1) + (1 — exp[z«;r]) + (1 — expfi/cr'])} 
— J d 2 K(f) WW (K)(l - exp[zKr]) / d 2 K(f) WW (K')(l - exp[i«V]) (79) 
so that the single-quark spectrum (}4*3*1 takes the form 

d^d~z = (^{/^W-)I(71P)-(71P--)| 2 
- J d 2 K(f) WW (K)({j*\p) - (7*|p - k)) 



(80) 



da D 



d 2 hd 2 pdz (27T) J 



J d 2 ^ w (^)(( 7 *|p) - ( 7 *|p - K ))| 2 . (81) 



Putting the inelastic and diffractive components of the FS quark spectrum together, we 
evidently find the FS parton density which exactly coincides with the IS parton density 
(|78j) such that p is indeed the transverse momentum of the FS sea quark. The interpre- 
tation of this finding is not trivial, though. 

Consider first the domain of p 2 < Q 2 < Q 2 A such that the nucleus is opaque for all 
color dipoles in the photon. Hereafter we assume that the saturation scale Q 2 A is so large 
that p 2 , Q 2 are in the pQCD domain and neglect the quark masses. In this regime the 
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nuclear counterparts of the crossing diagrams of figs. 3b,3d,3f can be neglected. Then, 
in the classification of [TT], diffraction will be dominated by the contribution from the 
Landau-Pomeranchuk diagram of fig. 3e with the result 



dq 



FS 



d 2 bd 2 p 



1 Q 2 da 



D 



D 



2 4ir 2 a em d 2 hd 2 p 

1 Q 2 

2 ' 4ir 2 a e 
N K 



■ J dz J d 2 K(f>ww{ K ) 



2 |(7*|P)| 2 



4vr 4 



12) 



Up to now we specified neither the wave function of the photon nor the spin nor the color 
representation of charged partons, only the last result in ()82j) makes an explicit use of 
the conventional spin-| partons. Remarkably, diffractive DIS measures the momentum 
distribution of quarks and antiquarks in the qq Fock state of the photon. We emphasize 
that this result, typical of the Landau-Pomeranchuk mechanism, is a completely generic 
one and would hold for any beam particle such that its coupling to hard colored partons 
is weak. In contrast to diffraction off free nucleons 01101 Ell- diffraction off opaque nuclei 
is dominated by the anti-collinear splitting of hard gluons into soft sea quarks, n 2 ^> p 2 . 
Precisely for this reason one finds the saturated FS quark density, because the nuclear 
dilution of the WW glue is compensated for by the expanding plateau. The result (f%2~j) 
has no counterpart in DIS off free nucleons because diffractive DIS off free nucleons is 
negligibly small even at HERA, r] D < 6-10 % [14J. 

The related analysis of the FS quark density for truly inelastic DIS in the same domain 
of P 2 < Q 2 < Q 2 A gives 



dq F s 
d 2 hd 2 v 



1 Q 2 f f 

2 ' 4n 2 a ' J Z 1 ^ K ^ WW ^ K "> ^ 7 *l p ~~ K ) 

Q 2 rQ 2 r 

— -= 4>ww{0) / d k d2|(7*|«) 



|2 



|2 



Nc Q2 fHO 2 n 2 ^ 



(83) 



It describes final states with color excitation of a nucleus, but as a function of the photon 
wave function and nuclear WW gluon distribution it is completely different from (JHj) for 
free nucleons. The ^-function simply indicates that the plateau for inelastic DIS extends 
up to p 2 < Q\. For Q 2 Q\ the inelastic plateau contributes little to the transverse 
momentum distribution of soft quarks, p 2 < Q 2 , but the inelastic plateau extends way 
beyond Q 2 and its integral contribution to the spectrum of FS quarks is exactly equal 
to that from diffractive DIS. Such a two-plateau structure of the FS quark spectrum is a 
new finding and has not been considered before. 

Now notice, that in the opacity regime the diffractive FS parton density coincides with 
the contribution oc | (7* |p) | 2 to the IS sea parton density from the spectator diagram of fig. 
3a, whereas the FS parton density for truly inelastic DIS coincides with the contribution 
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Figure 9: The two-plateau structure of the momentum distribution of FS quarks: (i) 
Q 2 ^ Q\a '■ The inelastic plateau is much broader than the diffractive one. (ii) Q 2 > Q 2 A 
: inelastic plateau is identical to diffractive one. 



to IS sea partons from the diagram of fig. 3c. In this regime the contribution from the 
crossing diagrams of figs. 3b, d is negigibly small. 

Our results (JH2~j) and (|H3*j) . especially nuclear broadening and unusually strong Q 2 
dependence of the FS/IS parton density from truly inelastic DIS, demonstrate clearly 
a distinction between diffractive and inelastic DIS. Our considerations can readily be 
extended to the spectrum of soft quarks, p 2 < Q 2 A) in hard photons, Q 2 > Q\. In 
this case the result (|82j) for diffractive DIS is retained, whereas in the numerator of the 
result (JHH1) for truly inelastic DIS one must substitute Q 2 — > Q\, so that in this case 
dqFs\o ~ dqFs\in and dqis ~ 2dqFs\o- The evolution of soft nuclear sea, p 2 < Q% is 
entirely driven by an anti-collinear splitting of the NSS-defined WW nuclear glue into the 
sea partons. 

The early discussion of the FS quark density in the saturation regime is due to Mueller 
Mueller focused on Q 2 ^> Q\ and discussed neither a distinction between diffractive 
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and truly inelastic DIS nor a Q 2 dependence and broadening (|83|) for truly inelastic DIS 
at Q 2 <Q\. 

A comparison with the IS nuclear parton densities which evolve from the NSS-defined 
WW nuclear glue shows an exact equality of the FS and IS parton densities. The plateau- 
like saturated nuclear quark density is suggestive of the Fermi statistics, but our principal 
point that for any projectile which interacts weakly with colored partons the saturated 
density measures the momentum distribution in the qq,gg, ... Fock state of the projectile 
disproves the Fermi-statistics interpretation. The spin and color multiplet of colored 
partons the photon couples to is completely irrelevant, what only counts is an opacity 
of heavy nuclei. The anti-collinear splitting of WW nuclear glue into soft sea partons 
is a noteworthy feature of the both diffractive DIS and IS sea parton distributions. The 
emergence of a saturated density of IS sea partons from the nuclear-diluted WW glue is due 
to the nuclear broadening of the plateau (jBHj) • Because the predominance of diffraction is 
a very special feature of DIS [2] , one must be careful with applying the IS parton densities 
to, for instance, nuclear collisions, in which diffraction wouldn't be of any significance. 



9 Signals of saturation in exclusive diffractive DIS 

A flat p 2 distribution of forward q, q jets in truly inelastic DIS in the saturation regime, 
see eq. (jHSJ), must be contrasted to the standard DGLAP spectrum, 

dq 2a< ? (p 2 )G(p 2 ) ^ 



d 2 p 37rp 2 

for the free nucleon target. 

In the diffractive DIS the signal of saturation is much more dramatic: a flat p 2 distri- 
bution of forward q,q jets in diffractive DIS in the saturation regime, see eq. (fS2|). must 
be contrasted to the spectrum 

da D G 2 (p 2 ) 

K w ( } 

in diffractive DIS off the free nucleon target |H1E3IHI|. Recall that it is p 2 which serves as a 
hard scale for diffractive dijets unless p 2 is so large that the pomeron splitting mechanism 
takes over [TTHI33]. 

In the general case one must compare the saturation scale Qa to the relevant hard 
scale for the specific diffractive process. For instance, in the exclusive diffractive DIS, i.e., 
the vector meson production, the hard scale for the proton target equals 



Q 2 «i(Q 2 + r4) (g6) 



and the transverse cross section has been predicted to behave as 

a,: X — — ( 87 ) 
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Figure 10: The modification of the {Q 2 + rriy) dependence of diffractive p production from 
a free nucleon target ( small saturation scale Q\) to a nucleus (large Q\) with the change 
of the saturation scale Q\. 



At Q 2 > Q\ the same would hold for nuclei too, but in the opposite case of Q 2 < Q\ the 
Q 2 - dependence is predicted to change to 



0~t OC 



G(x,Q 2 

Q\Q 2 



(88) 



The numerical results shown in Fig. 10 give an idea on how the (Q 2 + my) dependence 
of exclusive diffraction into the vector mesons changes from the free nucleon to nuclear 
target with the incresaing saturation scale. 

As shown in jHT], there is a duality correspondence between the (Q 2 +my) dependence 
of exclusive production of vector mesons and the mass spectrum of diffractive DIS into 
small-mass continuum states, i.e, at large values, (3 — > 1, of the diffractive Bjorken variable 



Q 2 

M 2 + Q 2 ' 



(89) 
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Specifically, for diffractive DIS off free nucleons the transverse diffractive structure func- 
tion has the large-/? behavior [33] 

F D oc (1-/3) 2 , (90) 
whereas in diffractive DIS off nuclei in the saturation regime we predict 

F D oc (1 - P) . (91) 

10 Jet-jet (de) correlation in DIS off nuclear targets: 
hard jets 

The diagonalization of the 2x2 matrix 04 derived in section 3 is a straightforward 
task, so that technically eqs. (fUjl and pH |) -(|3*0 jl allow a direct calculation of the jet- 
jet inclusive cross section in terms of the color dipole cross section a(r). The practical 
evaluation of the 6-fold Fourier transform is not a trivial task, though, and anlytical 
evaluations are called upon. Here we focus on production of forward hard jets with the 
momenta > Q\, where Qa is the saturation scale. Such a hard jets are produced from 
interactions with the target nucleus of small color dipoles in the incident photon such that 
diffractive nuclear attenuation effects can be neglected. In this case useful and intuitively 
transparent analytic results for the jet-jet (de) correlation can be obtained. 

Before proceeding to this approximation, it is convenient to introduce the average 
impact parameter 

b = i(b + + b' + + b_ + b'_), (92) 

(one should nor confuse b with the center of gravity of color dipoles where the impact 
parameters b ± and b'± must be weighted with z±) and 

s = b + - b^ , (93) 

for the variable conjugate to the decorrelation momentum, in terms of which 

b + - b'_ = s + r' , (94) 

b ^ b' + = s - r , (95) 

b - b'_ = s - r + r' . (96) 
We also reproduce here the matrix elements of 04 in this basis: 

<11|<t 4 |11> = tr(r) + o-(r') (97) 

(111(74188) = . 1 \a(s) - a(s + r') - a(s - r) + a(s - r + r')l (98) 
- 1 
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(88|a 4 |88) 



+ 



^_|[[ ff ( B ) + ff (8-r + rO] 
_J_[ a ( s + r ') + a ( s _ r )] 



1 



(T r 



a(r')] 



iVp — 1 

Hard jets correspond to |r|, |r'| <C |s|. Then the two eigenvalues are 



E 2 = (11|<t 4 |11) ~ 



and 



Ei = (88|ct 4 |88) « 2 



o-(s) = 2A c a(s) 



(99) 



(100) 



(101) 



where A c = N 2 /(N 2 — 1). Because of (|100|) only the second term, oc (ll|a 4 |88), must 
be kept in the Sylvester expansion (|3*T|). It is convenient to define the nuclear distortion 
factor 



exp 



exp 



-2 S 2T(b) 



(102) 



such that to the leading order in nuclear thickness, i.e., in the impulse approximation, 
Da{s,t,t') = 1. Then, the impulse approximation cross section per unit area in the 
impact parameter plane takes the form (we expose the relevant steps in some detail) 



daiA 



-1 



d 2 sd 2 rd 2 r' 



d 2 hdzd 2 p + d 2 p_ 2(2vr) 4 
x exp[-z(p + + p_)s + zp_(r' - r)]**(r')*(r) 



xT(b)y/N*-l(U\* 4 \88)D A (s,T,r f , b) 

= 2^)I T ( b ) / d 2 sd 2 rd 2 r' exp[-z(p + + p_)s + p_(r' - r)] 

xM/*( r ')\P(r) [cr(s) - cr(s + r') - cr(s - r) + a(s - r + r')] 

Here we employ the integral representation (jlj): 

[cr(s) - cr(s + r') - cr(s - r) + cr(s - r + r')] 
= a s a J d 2 Af(A) expfiAs] {1 - exp[zAr']} {1 - exp[-iAr]} , 

which leads to the Impulse Approximation result 

x exp[i(A — p+ — P-)s] exp[ip_(r' — r)] 
x\l/*(r') {1 - exp[iAr']} Hf(r) {1 - exp[-«Ar]} . 



(103) 



(104) 



(105) 
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The d 2 s integration entailes A = p + — p , the remaining Fourier transforms are straight- 
forward and give precisely (jl2J) times T(b): 

r ( b)— (106) 



d 2 bdzd 2 p + d 2 A d 2 bdzd?p + d 2 A 
The nuclear distortion factor takes a simple form 



cxp 



-^iT(b) 



1 ri 



D A (s,b) = j— =/ d/3exp 



J dp J d 2 K$(2/3A c z/ j4 (b),K)exp(zKs) (107) 



The introduction of this distortion factor into f)105j) is straightforward and gives our central 
result for the hard jet-jet inclusive cross section: 

d(Jin mfi \ f 1 m f i2 ^f n n\ /i \ a \ d<J N 



<Pb<HPr> <PA m Jo di3 1 A ~ ^ dz dC d 2 K ■ ( 108 ) 

(1 YJilZUj ]J-\-CL —A JO ./ (IZCL p + d re 

It has a probabilistic form of a convolution of the differential cross section on a free 
nucleon target with <3>(2/3A c z/,4(b), re). Here /3 has a meaning of the fraction of the nuclear 
thickness which the (qq) pair propagates in the color octet state. 
In the practical calculations one can use an explicit expansion 

/'»(2fe(HA) 'E ^^yS , (109) 
Jo t-^ V- W\rVA\b) 



3=0 



where j(J,x) is an incomplete Gamma-function. For heavy nuclei the dominant contri- 
bution comes from j ~ u A . Neglecting the small departure of A c from unity, for a heavy 
nucleus we can approximate 



f 1 dp$(2f3\ c u A (b), A) « $(\ c v A (b), A) 
Jo 



~ t{ ^ b) - A ^lwTQW (uo) 

and evaluate the gross features of a jet-jet decorrelation analytically. First, notice that 
/(re) which enters the free nucleon cross section (|12|) is a steep function of re compared 
to a broad distribution (jllOj) . Then, 

da m ^ Q\ da N 



d 2 bdzd 2 p + d 2 A v \{A 2 + Q 2 A ) 2 J dzd 2 p + d 2 n 
« l -T{b)e 2 f a em a s {p 2 + )[z 2 + (1 - z) 2 ]G{A 2 + Q*) 

x g| ^ 4 + (pD 2 fn r 

tt(A 2 + Q» )* (e 2 + plY ' 1 ' 
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where the integrated gluon structure function of the nucleon enters at the factorization 
scale ~ (A 2 + Q A ). Neglecting the scaling violations in G(A 2 + Q\), we can evaluate 
expectation value (A^_) subject to the condition A 2 < p+ as 




(H2) 



The differential out-of-plane momentum distribution can be evaluated as 



da 1 Q 



dA ± 2(Q 2 4 + Ai)3/2 



(113) 



which shows that the probability to observe the away jet at ~ disappears oc 1/Qa oc 

The probabilistic character of our result (|108|) can be understood as follows. Hard jets 
originate from small color dipoles. Their interaction with gluons of the target nucleus is 
suppressed by screening of color charges of the quark and antiquark in the color-singlet 
qq state which is manifest from the small cross section for a free nucleon target, see eqs. 
()12j) . ()13j) . The first inelastic interaction inside a nucleus converts the qq pair into the 
color-octet state. The color charges of the quark and antiquark are no longer screened 
and rescatterings of the quark and antiquark in the color field of intranuclear nucleons are 
uncorrelated. Consequently, the broadening of the momentum distribution with nuclear 
thickness follows the probabilistic picture. 



11 Disappearance of jet-jet correlation in DIS off nu- 
clear targets: minijets below the saturation scale 

In the above discusion of the single particle spectrum we discovered that the sea quarks 
evolve via the anticollinear, anti-DGLAP splitting of gluons into sea, when the transverse 
momentum of the parent gluons is larger than the momentum of the sea quarks jH] , which 
suggests strongly a complete azimuthal decorrelation of forward jets with the transverse 
momenta below the saturation scale, < Q\. Our numerical results reported in section 
6 suggest that for the realistic nuclei Q\ does not exceed several GeV 2 , hence this regime is 
a somewhat academic one. Still, let us assume that Qa is so large that jets with p^ < Q\ 
are measurable. 

Upon the slight generalization of (jl(J7J) the distortion factor admits a representation 
LU(s,r,r',b) = £d/3exp + (1 - /?)£ 2 ]T(b)} . (114) 

It is still too complex if one is after the exact eigenvalues £1,2, however, it simplifies 
greatly if one invokes the large- N c approximation, in which case 

S 2 = <ll|cr 4 |ll> = <r(r) + a(r') (115) 



31 



and 

E x = o-(s) + a(s + r' - r) . (116) 

Because in the large- N c approximation £2 = ( 1 1 1 cr 4 1 1 1 ) = a(r) + <7(r'), the first two terms 
in the Sylvester expansion (pTTf do vanish and only the last term will contribute to the 
jet-jet inclusive cross section. In the large-iV c approximation there is only one transition 
from the color-singlet to the color-octet state, once in the color octet state the qq only 
oscillates in the color space. Then the different exponentials in the resulting distortion 
factor 

D A (s, r, r', b) = J 1 d@ exp {-~/3[<t(s) + a(s + r' - r)]T(b)} 

x exp {--(1 - + a(r')]T(b)} (117) 

admit a simple interpretation: The last two exponential factors describe the intranuclear 
distortion of the incoming color-singlet (qq) dipole state, whereas the former two factors 
describe the distortion of the outgoing color-octet dipole states. 
A straightforward generalization of (jl(J7|) gives 

D A (s, r, r', b) = J d(3 j d 2 K^((l - P)v A (b), »i) exp(-m x r) 

x J d 2 K 2 ®((l - (3)v A (h), k 2 ) exp(iK 2 r) 

x J d 2 K 3 ®(f3v A (b), k 3 ) exp[irc 3 (s + r' - r)] 

x J d 2 K4$(f3v A (b), « 4 ) exp(zK 4 r) (118) 
and the jet-jet inclusive cross section takes the form 

d^bdtZ-dA = ^o, s a T(b) J^p J d^ ld ^ 2 d^ 3 d^f(.) 

X $((1 - (3)u A (b), Ki)$((l - P)V A Q>), K2) 

x <f>((3v A (b),K 3 )<l>(pv A (b),A- k 3 - k) 

X {^(-p_ + K 2 + #63) - ^(-P- + ^2 + ^3 + «)}* 

X {^(-p_ + Ki + K 3 ) - tf(-p_ + Ki + K 3 + «)} (119) 

It is uniquely calculable in terms of the NSS-defined WW glue of the nucleus. 

Now consider the limiting case of minijets with the transverse momentum below the 
saturation scale, |p_|,|A| < Q A . Notice, that ~ Q A , so that one can neglect p_ 
in the photon's wave functions and the decorrelation momentum A in the argument of 
$(/3v A (b), A — k 3 — k). The gross features of the product of the photon's wave functions 
which enters the integrand of after the averaging over Ki and n 2 can be well 

approximated by 

({# (-p_ + n 2 + k 3 ) - tf(-p_ + k 2 + k 3 + k)}* 

(~P- + Kl + K 3 ) - (-p- + Ki + K 3 + K)»K 1>2 

=}► |^( K3 )-^(k 3 + «)| 2 (120) 
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The principal point is that the minijet-minijet inclusive cross section depends on neither 
the minijet nor decorrelation momentum, which corroborates the anticipated complete 
disappearance of the azimuthal decorrelation of jets with the transverse momentum below 
the saturation scale. 

12 Summary and conclusions 

We formulated the theory of DIS off nuclear targets based on the consistent treatment 
of propagation of color dipoles in nuclear medium. What is viewed as attenuation in 
the laboratory frame can be interpreted as a fusion of partons from different nucleons 
of the ultrarelativistic nucleus. Diffractive attenuation of color single qq states gives 
a consistent definition of the WW unintegrated gluon structure function of the nucleus 
[TT] IH]. In these lectures we demonstrated how all other nuclear DIS observables - sea quark 
structure function and its decomposition into equally important genuine inelastic and 
diffractive components, exclusive diffraction off nuclei, the jet-jet inclusive cross section, - 
are uniquely calculable in terms of the NSS-defined nuclear WW glue. This property can 
be considered as a new factorization which connects DIS in the regimes of low and high 
density of partons. 

The two-plateau single parton (jet) inclusive cross section with the strong Q 2 depen- 
dence of the plateau for truly inelastic DIS has not been discussed before. A comparison 
with the initial state nuclear parton densities which evolve from the NSS-defined WW 
nuclear glue shows an exact equality of the FS and IS parton densities. The plateau-like 
saturated nuclear quark density is suggestive of the Fermi statistics, but our principal 
point that for any projectile which interacts weakly with colored partons the saturated 
density measures the momentum distribution in the qq,gg,... Fock state of the projec- 
tile disproves the Fermi-statistics interpretation. The spin and color multiplet of colored 
partons the photon couples to is completely irrelevant, what only counts is an opacity of 
heavy nuclei. 

The anti-collinear splitting of WW nuclear glue into soft sea partons is a noteworthy 
feature of the both diffractive DIS and IS sea parton distributions. The collective nuclear 
WW gluon field exhibits strong nuclear suppression, and the emergence of a saturated 
density of sea partons from the nuclear-diluted WW glue is due to the nuclear broadening 
of the plateau. Because the predominance of diffraction is a very special feature of DIS 
j2], one must be careful with applying the IS parton densities to, for instance, nuclear 
collisions, in which diffraction wouldn't be of any significance. The anti-collinear splitting 
of WW nuclear glue is a clearcut evidence for inapplicability of the DGLAP evolution to 
nuclear structure functions unless Q 2 ^> Q\a- 

Although in our derivation we focused on DIS, all the results for hard single-jet and jet- 
jet inclusive cross sectinos are fully applicable to production of jets in the beam fragmen- 
tation region of meson-nucleus collisions and can readily be extended to nucleon-nucleus 
collisions. Indeed, as argued in [TT], the final state interaction between the final state 
quark and antiquark can be neglected and plane-wave approximation becomes applicable 
as soon as the invariant mass of the dijet exceeds a typical mass scale of prominent meson 
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resonances. 

Of particular interest are the results on the azimuthal decorrelation of hard jets, in 
particular, disappearance of azimuthal jet-jet correlations of minijets with momenta below 
the saturation scale. For the average DIS on heavy nuclei we estimate Ql A ~ 2 GeV 2 , but 
for the central DIS at a small impact parameter the saturation scale can be several times 
larger. For instance, from perihperal DIS to central DIS on a heavy nucleus like Pt, i>sa 
rises form u 8 a = 1 to u 8 a ~ 13, so that according to (|113|) a probability to find the away 
jet decreases by the large factor ~ 3.5 from peripheral to central DIS off Pt target. 

The early experience with application of color dipole formalism to hard hadron-nucleus 
interactions |21. suggests that our analysis can be readily generalized to mid-rapidity jets. 
One only has to choose an appropriate system of dipoles, for instance, the open heavy 
flavor production can be treated in terms of the intranuclear propagation of the gluon- 
quark-antiquark system in the overall color-singlet state. For this reason, we expect similar 
strong decorrelation of mid-rapidity jets in hadron-nucleus and nucleus-nucleus collisions. 
To this end, recently the STAR collaboration reported a disappearance of back-to-back 
high p± hadron correlation in central gold-gold collisions at RHIC flj . Nuclear enhance- 
ment of the azimuthal decorrelation of the trigger and away jets discussed in sections 
10 and 11 may contribute substantially to the STAR effect. Here we emphasize that 
we discuss the distortions of the produced jet-jet inclusive spectrum due to interactions 
with the nucleons of the target, a practical consideration of azimuthal decorrelations in 
central heavy ion collisions must include rescatterings of parent high-p^ partons on the 
abundantly produced secondary hadrons. 

This work has been partly supported by the INTAS grants 97-30494 & 00-00366 and 
the DFG grant 436RUS17/89/02. 
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